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ABSTRACT 

Context. CoRoT, the pioneer space-based transit search, steadily provides thousands of high-precision light curves with continuous time sampling 
over periods of up to 5 months. The transits of a planet perturbed by an additional object are not strictly periodic. By studying the transit timing 
variations (TTVs), additional objects can be detected in the system. 

Aims. A transit timing analysis of CoRoT-lb is carried out to constrain the existence of additional planets in the system. 

Methods. We used data obtained by an improved version of the CoRoT data pipeline (version 2.0). Individual transits were fitted to determine the 
mid-transit times, and we analyzed the derived O - C diagram. N-body integrations were used to place limits on secondary planets. 
Results. No periodic timing variations with a period shorter than the observational window (55 days) are found. The presence of an Earth-mass 
Trojan is not likely. A planet of mass greater than ~ 1 Earth mass can be ruled out by the present data if the object is in a 2:1 (exterior) mean 
motion resonance with CoRoT-lb. Considering initially circular orbits: (i) super-Earths (less than 10 Earth-masses) are excluded for periods less 
than about 3.5 days, (ii) Saturn-like planets can be ruled out for periods less than about 5 days, (iii) Jupiter-like planets should have a minimum 
orbital period of about 6.5 days. 
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1. Introduction 

As a consequence of the gravitational perturbations, the mid-times of 
consecutive transits deviate from a linear ephemeris in a transiting ex- 
oplanet system (transit timing variation, hereafter TTV). Depending on 
the mass and the orbital configuration of the perturbing planet, this de- 
viation can have amplitudes from a few seconds to days (e.g. Steffen, 
2006). Moreover, the duration, shape, and depth of the transits can also 
change. In extreme cases the transits can disappear and then reappear 
(Schneider 1994, 2004). Both the theoretical aspects and the observ- 
able effects have been studied in e.g. Miralda-Escude (2002), Borkovits 
et al. (2003), Agol et al. (2005), Holman & Murray (2005), Ford & 
Holman (2007), Simon et al. (2007), Heyl & Gladman (2007), Steffen 
et al. (2007), Nesvorny & Morbidelli (2008), Pal & Kocsis (2008) and 
Kipping (2009). Several transiting exoplanets were subject to this kind 
of analysis (Steffen & Agol 2005, Agol & Steffen 2007, Miller-Ricci 
et al. 2007, Alonso et al. 2008, Hrudkova et al. 2008, Miller-Ricci et 
al. 2008ab, Diaz et al. 2008, Coughlin et al. 2008, Rabus et al. 2009, 
Stringfellow et al. 2009). 

In addition, stellar spots may affect the transit shape and, because of 
this, we have some difficulty in determining of the midtime of the tran- 
sit. This effect may cause spurious periodic terms in the O — C diagram 
of exoplanets (Alonso et al. 2008, Pont et al. 2007). 

Here we report the TTV analysis of CoRoT-lb based on data ob- 
tained by an improved version of the CoRoT data pipeline. In this sys- 
tem a low-density planet (mass: l.03M Jllp , radius: l.49Rj up , average 
density: 0.38 gcm~ 3 , semi-major axis: 5A6R e , orbital period: 1.5089557 
days) orbits a G2 main-sequence star (Barge et al. 2008). Thirty-six 
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transits were observed by CoRoT, 20 of them in 512 seconds and 16 
with the 32 seconds sampling rate mode. In total more than 68 000 data 
points were collected during 55 consecutive days (Barge et al. 2008). 
The operation of the satellite is described in detail in the pre-launch 
book, and the reader can find useful information about CoRoT in Baglin 
et al. (2007), Boisnard et al. (2006), and Auvergne et al. (2009). 

2. Methods of mid-transit point determination: 
effect of the sampling rate on the precision 

If one uses the CoRoT data for TTV analysis, the main limiting factor 
arises from the sampling rate. The typical length of the ingress/egress 
phase of a hot Jupiter is on the order of 10-20 minutes. (In the particular 
case of CoRoT-lb, the ingress/egress time is 9.8 minutes.) CoRoT tar- 
gets are observed with 512 or 32 second sampling rates (the so-called 
undersampled/oversampled modes, see Fig 1). Concerning a typical 3 
hour transit, one can easily conclude that a transit observation consists 
of only (3 x 3600seconds)/512seconds » 21 data points. In the over- 
sampled mode we have typically over 300 data points per transit. The 
small number of data points in the undersampled mode may not be bal- 
anced by the very good photometric precision of CoRoT (which is about 
0.1% for a 13 magnitude star in white light for a 512 second exposure, 
see Auvergne et al. 2009), therefore we chose to investigate this issue. 

A second factor arises from the orbit of the satellite. The satellite 
periodically crosses the so-called South Atlantic Anomaly (SAA) re- 
gion, which causes bad/uncertain data points and long data gaps (typi- 
cally 10 minutes). This is significant only when the SAA-crossing oc- 
curs during the ingress or the egress phase. Therefore the following test 
was carried out. Using the exoplanet light curve model of Mandel & 
Agol (2002) and the parameters of the system (Barge et al., 2008), we 
simulated the light curve of CoRoT-lb. Then this curve was re-sampled 
to the same time-points as CoRoT observations. We added a Gaussian- 
like random noise to the points. The standard deviation of the noise term 
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Table 1. Mid-transit times (HJD) and errors (days) of CoRoT-lb. 
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was chosen in such a way that we had the same signal-to-noise ratio as 
given in Barge et al. (2008) for the CoRoT-1 light curve. A constant 
orbital period was assumed. 

Then we determined the mid-transit times in this simulated light 
curve, fitting each individual transit separately. Again, we use the 
Mandel & Agol (2002) model combined with the Amoeba algorithm 
(Press et al. 1992) to find the optimum fit. We assumed that the planet- 
to-stellar radius ratio and the limb-darkening coefficients are known, so 
they were fixed. The adjustable parameters are the mid-transit point, the 
inclination, and the a/R s ratio (a: semi-major axis, R s : stellar radius). 

On average, the fits of the individual transits yield only a difference 
of 9 seconds between the real and the determined midtransit points in 
the undersampled mode, when there are data points in the ingress/egress 
part of the light curve. When the ingress or egress part is missing in 
the undersampled mode, the errors can be as large as 20-60 seconds, 
depending on the distribution of points during the transit. If both the 
ingress and egress parts are missing, the errors are 60-120 seconds, 
sometimes even more. 

These optimistic error bars should be increased due to at least two 
different effects. First, we do not know a priori the exact value of the 
orbital period which leads to small uncertainties in the calculation of 
the phase. We estimate that for CoRoT-lb this is negligible. Second, 
the stellar activity is not included. However, the detailed discussion of 
these two effects goes beyond the purpose of the present investigation. 

Since we use a constant period and assume e=0 during the sim- 
ulation, we expect a linear O-C curve with some scatter. To better 
characterize this scatter we calculate the standard deviation of the sam- 
ple. We find that the mean lcr scatter of the resulting overall O-C 
diagram of this simulated light curve with a constant period is 22 sec- 
onds. But it is 27 seconds for the undersampled part and 16 seconds for 
the oversampled part. 



3. TTV analysis of CoRoT-1 b 

We used the N2 level data points (Auvergne et al. 2009) processed by 
the 2.0 version of the pipeline (not yet realeased data for the public). 
The resulting light curve was manually checked: a few data points were 
noted by the pipeline to be affected by cosmic ray events in spite of it 
having no problems - we restored these data points. In addition, several 
outliers were removed by hand. Then we performed a TTV analysis by 
fitting all transits using the method described in the previous section. 
Transit No. 30 is excluded from this investigation because it is strongly 
affected by noise. Table 1 gives the midtransit times and their errors. 

The overall O-C diagram and its Fourier-transform are given in 
Fig. 2. This diagram is built using the observed light curve. It is promi- 



nent that after switching on the 32 sec sampling rate mode (after the 
20th transit), the lcr scatter of the O-C diagram is reduced to only 18 
seconds, compared to the lcr scatter of the 25 seconds observed in the 
undersampled mode. The lcr scatter of the whole O-C diagram is 22 
seconds. All these scatter values are very close to the value we would 
expect in the case of a constant orbital period (see previous section). 

No clear periodicity or trend can be identified in this diagram. We 
calculate the Fourier-spectrum of the O-C curve by the Period04 soft- 
ware (Lenz & Breger 2005) to search for any non-obvious periodicity. 
The power spectrum shows few peaks, but none of them is above the 
noise level. The highest peak has only S/N x 2, so it is not regarded as 
a real signal. 

We conclude that there are not any periodic TTVs in CoRoT-lb 
with a period less than the observational window (55 days) and an am- 
plitude larger than about 1 minute (=3<r detection level). We also note 
that there is no significant change in the inclination and the a/R star ratio 
during this interval. Bean (2009) presents results about the TTV analy- 
sis of CoRoT-lb based on the public data processed by an earlier version 
of the pipeline. His O-C diagram shows a larger ltr-scatter (36 seconds 
compared to our 22 seconds). 

4. Limits on secondary planets 

4. 1 . Examples of simulated TTVs 

We show in Table 2 which amplitudes and periods can typically be 
expected in the O-C diagram of the CoRoT-lb system based on 
some dynamical simulations. We include an Earth-mass planet at the L 4 
Lagrangian point, an Earth in the 2: 1 exterior mean motion resonance, 
a nearby Neptune-like planet, or an outer eccentric Neptune. The stel- 
lar mass, the mass of CoRoT-lb and its orbital elements are fixed to 
the values given in Barge et al. (2008). We use the Mercury software 
(Chambers 1999), with the Burlish-Stoer algorithm (accuracy parame- 
ter 6 = 10~ 16 ). As one can see in Table 2, O - C variations on the order 
of 60-150 seconds might occur on a short time scale (typically 10-150 
orbital cycles of the transiting planet). 

Ca.se 1: an Earth-mass Trojan planet librating with an amplitude 
of 20° around the L 4 Lagrangian point would have an amplitude of 60 
seconds in the O-C diagram with a period of about 10 orbital cycles 
of the transiting planet (about 15 days, see Table 2). The amplitude is 
close to our detection limit. There is a peak in the Fourier- spectrum of 
the O-C diagram at the corresponding frequency with an amplitude of 
about 1 1 seconds. However, the peak is not significant (S/N =1.3 only). 
Therefore, an additional planet with similar parameters is not likely. 

Ca.se 2: An Earth-mass planet initially on a circular orbit and in 
2:1 mean motion resonance with CoRoT-lb would have an amplitude 
of about 100 seconds in the O-C diagram with a period of about 150 
transits (approximately 225 days, see Table 2 and Fig. 3). If the CoRoT 
observational window was around the maximum or the minimum of the 
O-C curve (see Fig. 3) then we would have no chance to discover this 
possible planet because the amplitude is on the order of the scatter. If the 
observational window matched the steepest part of the O - C diagram, 
we would observe a linear O-C curve that could be interpreted as a 
wrong period value. This gives a hint: if there are no observed period 
variations in a short observational window, this does not mean that we 
can give an upper limit for a hypothetical perturber object. It might be 
the case that we are on a linear part of the O-C curve. The observational 
window should be long enough to exclude similar cases. 

Ca.se.s- 3 and 4: simulations show that an outer 30 Earth-mass planet, 
close to CoRoT-lb (P = 2.772118632 days and e = 0.05) or eccentric 
CP = 4.2679123 days and e = 0.25), cause O-C variations of about 150 
seconds, within approximately 15 and 30 orbital revolutions of the tran- 
siting planet, respectively (see Table 2). This is much greater than our 
detection limit, so outer planets in the CoRoT-lb system with similar 
orbital parameters can be excluded. 

4.2. Detailed analysis 

Using N-body integrations, we computed the maximum mass of a hy- 
pothetical perturbing planet, with given initial orbital periods and ec- 
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Fig. 1. Top: A typical transit observation of CoRoT-lb by CoRoT in the 512 sec sampling rate mode. Bottom: A typical transit 
observation of CoRoT-lb by CoRoT in the 32 sec sampling rate mode. Abcissa is phase, ordinate is normalized intensity. The solid 
lines show the fit. 
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Fig. 2. Top: The O-C diagram of CoRoT-lb obtained by the light curve fit of the individual transits. Bottom: The power spectrum 
of the Fourier-transform of the O - C diagram. Frequency is in cycles/day and amplitude is in seconds. 



centricities, leading to TTVs compatible with the data. To calculate the 
transit times, we used a bracketing routine from Agol et al. (2005). The 
orbits of CoRoT-lb and an additional planet were computed over the 
timespan of the observations, using a Burlish-Stoer integrator with an 
accuracy parameter 6 = 10~ 16 . The equations of motion were integrated 
in a cartesian reference frame centered on the barycenter of the system. 
The transit times are subtracted from the data to give the O-C residuals 
and^ 2 . 

The masses of the central star and CoRoT-lb are respectively fixed 
at 0.95 solar masses and 1.03 Jupiter masses (Barge et al. 2008). The 
orbit of CoRoT-lb is initially circular with an orbital period P = 
1.5089557 d (Barge et al. 2008) and a true longitude 9 = deg. With 
these parameters and without any perturbation due to an additional 
body, the first transit occurs at T(HJD) = 2454138.327840, and the 
residuals given by the numerical integration are at their minimum (i.e. 
the same as the ones from the best constant period fit, see Sect. 3) with 
the following values: zero mean, standard deviation <x mi „ = 21.62 sec- 
onds,^ 2 = 24.55. 

The perturbing planet is assumed to lie on the same plane as 
CoRoT-lb. For given initial orbital parameters, we increase the mass 
of the test planet, starting at 0.1 Earth masses, and calculate the stan- 
dard deviation of the O-C residuals. We store the mass value for which 
this rms exceeds the observed scatter cr min . In this way we determine 
the maximum planet's mass allowed. The results are shown in Fig. 4 
(respectively Fig. 5), which shows the maximum mass for a perturb- 
ing object as a function of its initial orbital period (resp. initial orbital 
period and eccentricity). In Fig. 4, the mass of the secondary planet 
has been varied between 0.1 and 100 Earth masses (100 values on a 
log scale), and its initial orbital period between 2.8 and 7.6 days (with a 
step of 0.0015 days). For any given orbital period, eccentricity, and mass 
value, the TTV-signal is computed over the range of possible initial true 
anomaly and longitude of pericenter values to minimize the resulting 



residuals. In Fig. 5, the perturbing planet is initially at its apocenter 
(fixed at 180 degrees from CoRoT-lb), and the following grid of pa- 
rameters has been used: (i) masses between 0.1 and 100 Earth masses 
(100 values on a log scale), (ii) orbital periods between 2.8 and 7.6 days 
(with a step of 0.001333 days), (iii) eccentricities between and 0.25 
(100 values on a log scale). 

From Fig. 4, Saturn-like planets can be ruled out for periods less 
than about 5 days if e=0 (respectively 6 days if e=0.2). As shown in the 
figures, perturbing planets with eccentric orbits obviously cause larger 
TTVs, hence have lower mass limits. Super-Earths are defined as plan- 
ets with 1-10 Earth masses (Valencia et al. 2007). Depending on the ini- 
tial eccentricity, such planets with orbital periods less than 3.4-4.1 days 
can be excluded. Planets with masses greater than 0.3-1.0 Earth masses 
can be ruled out by the data if they are in the 2: 1 (exterior) mean motion 
resonance with CoRoT-lb. The data do not allow strongly constraining 
the mass of perturbing planets near higher order resonances. Finally, we 
estimate the minimum orbital period for an outer Jupiter-mass planet. 
From Holman and Murray (2005), 

u - HA ^max I ^perturber \ 3 

™perturber - "TT" '"star p — I (I - e p crtiiiW Uj 

40 -(transiting \ 'transiting ) 

When assuming a circular orbit and At max = 3tr m ,„, this yields a mini- 
mum orbital period of 2.0 days. Otherwise, we would see its effect in 
the O - C diagram. This lower limit is in good agreement with the nu- 
merical simulations (see Figs. 4 and 5) 

5. Summary 

Our work shows that CoRoT allows study of the short time scale (30 
days for the Short Run fields, 150 days for the Long Run fields) tran- 
sit timing variations whose l<x detection limit depends on the sampling 
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Table 2. Amplitudes and periods of O — C variations in CoRoT- 
lb system. 
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~ 15 
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outer planet initially 
with e=0.25 and 
P=4.2679123 days 


150 


~ 30 



rate, and it is 22 seconds for CoRoT-lb. The comparison of the O — C di- 
agram of CoRoT-lb with numerical integrations leads to the following 
results: (i) an Earth-mass planet at the L4 point is not likely. If exist- 
ing, its detectability would be close to the 3<x detection limit, (ii) an 
outer Earth-mass planet in 2:1 resonance with CoRoT-lb can be re- 
jected, given our data set. However, a longer observational window is 
required to fully assess the presence of such a planet, (iii) super-Earths 
are excluded for periods less than about 3.5 days, (iv) Saturn-like plan- 
ets are ruled out for periods less than about 5 days. 

Bean (2009) finds that there is no additional planet in the system 
with 4 Earth-mass or greater on an orbit with 2:1 mean motion res- 
onance. Using an improved version of the C0R0T data pipeline, we 
confirm his result. 

We also showed that TTV analyses of C0R0T data are promising 
for detecting additional objects in transiting systems observed by the 
satellite. 
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Fig. 3. The simulated O-C diagram of CoRoT-lb if the transiting planet is perturbed by an Earth-mass planet initially on a circular 
orbit in 2:1 mean motion resonance. 
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Fig. 4. Maximum allowed mass of a hypothetical perturbing object as a function of its orbital period for excentricities e=0 and 0.2. 
The 2: 1 mean motion resonance is indicated. 
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Fig. 5. Upper mass limits for a hypothetical second object in the CoRoT-lb system as a function of the perturber's orbital period 
and eccentricity. 



